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SUMMARY 

Analysis of high -altitude -rocket electron and ion density measure- 
ments suggests a simple model of the quiet ionosphere. Near h m0JC F 2 
the profile is given by the "a-Chapman" function, that is, the Chapman 
electron density equation with sec x - 1 and the scale height H a con- 
stant. A method is given for determining h max F 2 , N max F 2 , and H from 
N(h)data obtained from ionograms. Well above the peak, the profile is 
taken to have a constant exponential slope of 200 km during the day and 
150 km at night. If simultaneous nearby measurements of the total 
electron content are available, a more accurate slope may be computed 
to provide the profile up to about 1000 km altitude. 
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A MODEL OF THE QUIET IONOSPHERE* 


by 

J, Carl Seddon 
Goddard Space Flight Center 


INTRODUCTION 

Electron density profiles above the f 2 maximum obtained with rockets and satellites are sum- 
marized by Wright in Reference 1. He found that he could obtain an approximate agreement with these 
results if he used the Chapman electron density equation with scale height H considered to be constant 
at 100 km and with sec x = 1. He computed the ratio of electron density above the maximum to that 
below; but these values tended to be a little higher than the measured values, which were obtained 
mainly from orbiting satellites and military rockets where accuracy was "not of the best." Sub- 
sequently Berning (Reference 2) obtained a sunrise profile to 1500 kilometers which, when normalized, 
did not agree well with previous profiles. His results indicated the existence of a considerable gra- 
dient in the electron scale height; however, this profile was — due to unforeseen circumstances 
obtained at ground sunrise when conditions could be changing at a rapid rate. A little later, Bowles 
(Reference 3) published a daytime profile obtained by the incoherent backscatter technique, which 
indicated an approximate Chapman distribution with H ~90 km. His night profile indicates H ^75 km, 
although the densities close to the maximum are somewhat larger than a Chapman distribution with 
this scale height. 

Nisbet and Bowhill (Reference 4) published a series of profiles obtained with military rockets 
under difficult scientific conditions. They attempted to compare their normalized results with a 
Chapman distribution by utilization of a variable neutral scale height given by Kallman in Reference 5; 
the agreement was not particularly good. Berning t revealed profiles showing a constant scale height 
Hwell above h n>x F 2 of about 100 km in the daytime and 72 km during the evening. 

Hanson and McKibbin (Reference 6) reported ion density measurements which gave an h value of 
75 km in the evening for altitudes well above the maximum density. Pineo et al.t reported on daytime 
results obtained by incoherent backscatter techniques, which showed a scale height gradient similar to 
the sunrise results of Berning. More recently, Jackson and Bauer (Reference 7) obtained a daytime 
profile with H also given as 100 km. 


•Presented at the URSI Ionosphere World-Wide Soundings Symposium, Nice. France, December 11-16, 1961. 

fBerning, W. W., ’Experimental Measurement of Electron Temperature and Electron and Ion Density by a Direct Measurement and a Propaga- 
tion Technique’.’ Paper presented at Union Radio Scientifique Internationale Meeting, Washington, May 1961. 

IPineo, V. C., et al., ’Expected Studies of F-Region Ending in Coherent Back Scattering Technique at Frequencies Around 400 Me. Paper 
presented at Union Radio Scientifique Internationale Meeting, Washington, May 1961. 
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This paper discusses the electron and ion density results obtained with rockets under quiet iono- 
spheric conditions and reasonably favorable scientific conditions. A simple model of the quiet iono- 
sphere, which can be expressed approximately in analytical form, is obtained. This model is used to 
develop a means of obtaining from ionograms an approximate electron density profile and total elec- 
tron content. It is also shown how such data used in conjunction with total electron content measure- 
ments make possible the determination of the electron density profile above \ ax ¥ 2 . 


DISCUSSION 

If the results ofBerning* , Hanson and McKibbin (Reference 6), and Jackson and Bauer (Reference 
7) are examined near the peak of the F 2 region, it is found that the measured electron densities do not 
follow a Chapman function with a scale height as high as 100 km. However, these results are all in 
agreement that, well above the peak altitude, a constant exponential slope exists which— on the assump- 
tion that the ion and electron temperatures are equal— corresponds to an h value of about 100 km in the 
daytime and 75 km in the evening. 

Yonezawa (References 8 and 9) showed that under the influence of vertical diffusion and nighttime 
attachment any N(h) curve will at night tend to take the form of a Chapman layer of constant scale 
height H and sec x = 1, referred to as an M a-Chapman f1 layer. Long (Reference 10) demonstrated that, 
if the N(h) analysis of the ionogram includes an allowance for the underlying ionization that is usually 
neglected, the profiles for all latitudes follow an a-Chapman variation to first order. The present 
paper shows that, for the limited daytime rocket data available, this also holds for quiet daytime 
conditions. 

Figure 1 gives the Jackson- Bauer results plotted near the F 2 peak. The circled points represent 
the a.- Chapman function with a scale height of 57 km, and the numbers show the number of scale 
heights above or below the maximum. There is indeed a very close agreement with this Chapman 
function over the range from - 1 scale height up to about + 1.6 scale heights with h max = 287 km. At 
380 km a break occurs in the profile slope until, at an altitude of about 400 km, a constant asymptotic 
slope is reached. This disagrees with the result given by Jackson and Bauer that hydrostatic equi- 
librium begins at 350 km. Figure 2 shows a portion of Berning f s unpublished quiet day results, with 
equilibrium also beginning at 400 km. His results follow a Chapman distribution very accurately from 
about - 1 to + 1.5 scale heights with H = 60 km and h max = 280 km. Figure 3 gives the nighttime ion 
density results of Hanson and McKibbin, which also show equilibrium beginning at 400 km. Their 
results follow a Chapman distribution from below -1.2 to +1.3 scale heights with H = 43 km and 
h ma% = 309.5 km. Berning f s unpublished night flight at Wallops Island, Va., at 2143 EST July 13, 
1960, also follows the Chapman function very accurately from - 2.0 to +1.8 scale heights with H = 63 
km and h max = 357 km. 

It is difficult to state precisely at what altitude the neutral particle scale height has a value equal 
to H, but it would seem reasonable to suppose that it is near h max F 2 . Kallmann f s (Reference 11) 


*See footnote, page l. 
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NEUTRAL SCALE HEIGHT (KALLMANN) (km) 






average daytime scale heights are shown on the 
right-hand side of Figure 1; the agreement is 
even better when compared with Berning’s day- 
time flight. Kallmann’s nighttime average value 
agrees reasonably well with Berning’s summer 
night firing, but rather poorly with Hanson and 
McKibbin’s November night firing. 

These data thus indicate that under quiet 
conditions the F 2 is a-Chapman from -2.0 to 
+ 1.5 scale heights at night and from -1.0 to 
+ 1.5 in the daytime. Above 1.5 scale heights 
there is a transition region of about 25 kilo- 
meters after which diffusive equilibrium exists, 
with a constant exponential slope of about 150 km 
in the evening and 200 km in the daytime. Data 
late at night are not available, and it is likely 
that before sunrise the slope may decrease to 
lower values; a method will be suggested to 
check on this. 

MODEL FOR THE QUIET 
MIDDLE-LATITUDE IONOSPHERE 

Since the ionosphere decreases the track- 
ing accuracy of systems using radio waves, it 
would be desirable to express the density pro- 
file in terms of approximate analytical functions 
to simplify calculations. Therefore the ioniza- 
tion in the E- region and above was reviewed, 
and a reasonable approximation in the daytime 
below -1.0 scale height may be had by drawing 



POSITIVE ION DENSITY [(ions/cm 3 ) x 1Q- 5 | 


Figure 3 — Chapman distribution fitted to Hanson and 
McKibbin ion density profile, (Nov. 9, I960, 2044 EST 
Wallops Island, Va.) 


a straight line on semilog paper to N max (E)at an altitude of 100 km. While the actual profiles indicate 
variations around this assumed profile, the electron contents are about equal. Figure 4 shows the ap- 


proximate model for the quiet middle- latitude ionosphere. 


Table 1 presents an approximate representation of the profiles in analytical form, which may be 
useful in refraction calculations; S is obtained from a knowledge of the initial and final values. Also 
given are formulas for the electron content. It is assumed here that the nighttime region is entirely 
Chapman below h max . If appreciable sporadic-E exists, it may be assumed to have a plasma frequency 
equal to the lowest F- reflection frequency plus one-half the gyrofrequency and an average thickness of 
1 km (Seddon, Reference 12). 
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Table 2 gives the electron content calcu- 
lated from the formulas of Table 1 as compared 
with the values obtained from numerical inte- 
gration of rocket data, with the assistance in 
some cases of E-region data obtained from the 
P'-f analysis performed by the National Bureau 
of Standards. 

Comparison of the rocket results and the 
results obtained by the National Bureau of 
Standards using ionograms from a nearby iono- 
sonde showed that the agreement was generally 
good but that considerable disagreement existed 
at times for H and/or h max . The National Bu- 
reau of Standards technique involves a parabolic 
assumption near h max that should give fairly 
good results if the ionogram can be read to 
frequencies very close to f Q F 2 - Below about 
-1/2 scale height, however, the parabolic 



middle-latitude ionosphere 


Table 1 


Analytical Representation of the Model 


Night N 

N_ (F a ) Ch(z) 

h - h 

z < 1.5 , z = 

H 

h - h „„x->- 5H 

0. 7N (F,) e 150 

max V 2/ 

h>h m „* +1.5H 

Day N 

h- 1 00 

N a *(E)e s 

N ma> (E)<N<0.7N m „ (F 2 ) , 

100 < h < h - H 

max 

N max (F 2 ) Ch(z) 

-1 < z < 1.5 

0-™ ma ,(F 2 )e 

h > h + 1.5H 

— - max 

Night 

N b 

1.312 HN (F,) 

max \ 2/ 

N a 

1. 30 HN (FJ| +150X0.7N (F„\ 

max V 2/ m a x V 2/ 

N t 

2.61 HN (FA + 150 x 0. 7N (F ) 

max\ 2 / max \ 2 / 

Day 

N b 

S l~0. 7N (FA -N (E)l + 0. 88 HN (FA 

N„ 

1. 30 HN max (F 2 ) + 200 x 0.7N max (F 2 ) 

N t 

s[0.7N raa)t (F 2 ) -N max (E)] + 2.18HN max (F 2 ) + 200 x 0. 7N max (F a ) 


5 





Summary of Results 



0 £ o M hS JS, 

r * r, 8 8 . 5 . 


S. 8.1 8.X 8. 8. S. 

O OjC 0 t5 O _0 


O IT T'fl' it 

o ^ 

v\ ■— ( -h O O 0 s - 

H N N N N O 

lift!! 
k\ -<r kr. « 


* I * 


>. > >> 

2 . 3 £ 



6 


JCbapumn function above N{h) reduced data. 



approximation becomes rapidly poorer; and examination of their data frequently shows few or no 
points higher than —1/2 scale height. In addition, the method depends on an accurate value of f Q F 2 . 

The same N(h) data were used to obtain the best fit possible to a Chapman distribution; and these re- 
sults generally gave better agreement with the rocket results including N mflX F 2 .which is obtained with- 
out the use of f Q F 2 . 

EXPERIMENTAL METHOD 

While various methods can be devised to fit theN(h) values obtained, the method used in this 
paper can be performed on a desk computer or an electronic computer. Both computers actually 
were used, the electronic computer being the relatively slow LG-30 digital computer. 

All the values obtained from theN(h) reduction that were above 1 scale height below the esti- 
mated h max in the daytime, and above about 1.5 scale heights below h max at night, were used. The 
highest altitude was called h v the next h 2 , etc. If the points fall on a Chapman curve, the relation 


must be constant for each value of i, where z is the solution of Ch(z) with sec x - 1; that is, 

N = N ma, exp | (1 - Z - e- z ) . (2) 

The values of z can be found from the relation 


Ni 

Ch( z . ) = • (3) 

max 

The computer is given a value of N max a little less than the value computed from the observed 
f F 2 . The h values are obtained from the N(h) reduction, along with the corresponding value of n. 

The result obtained is a series of values H 12 , h 13 , etc. which will increase steadily in value if the 
N iTiax used is too small and if the N(h) values are, in fact, truly Chapman. Because of the fact that 
there is some variation from a true Chapman function, the variation will not be entirely smooth. The 
mean value H is computed, and also the average deviation from the mean. 

A slightly larger value of N max is then given to the computer, and the same computations are made 
again. If the N(h) points were exactly on a Chapman curve, a value of N max ultimately would be found 
where the average deviation would be zero. Because of irregular variations from true Chapman, the 
average deviation decreases rapidly and nearly linearly to a value dependent on the size of such vari- 
ations and then increases again, as illustrated by Figure 5. The altitude h max F 2 is the average of the 
values obtained for h for the correct N , where 

max. max' 

i 

h ,., . = h i ” S z i • (4) 
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RESULTS 
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The values obtained for H, h max , and N max 
are shown in Table 2 in the columns labeled 
"Chapman," and those obtained by the NBS 
using the parabolic approximation are shown 
in the other columns. The table is arranged 
to facilitate comparison with the rocket re- 
sults. The Wallops Island ionosonde is only a 
few miles from the rocket launcher. The Fort 
Belvoir ionosonde is approximately 160 km to 
the west, or about 8 minutes earlier in solar 
time. 
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Figure 5— Determination of N max (F 2 ) by Chapman-fit to 
ionogram N(h) reduction. (April 27, 1961, 1500 EST Wal- 
lops Island, Va«) 


The values for both H and N max obtained by 
the Chapman technique generally agree better 
with the rocket results than those obtained 
from the parabolic approximation except for 
the one case of the winter night firing. It is 

quite likely that the ion trap measurements are not as accurate as the electron density measurements, 
since the N(h) reduction disagrees by a variable amount up to nearly 10 percent, which can introduce 
a considerable change in the value of H. 


In all the other cases examined, it was found that the N(h) reductions were in good agreement 
with the rocket data except that, above 1 scale height below the maximum, the values of N were always 
slightly lower. However, this slight decrease tends to make both H and h max F 2 too large and N max F 2 
slightly too small. The result causes N b to be too large in the daytime. A better result is obtained 
by integrating the N(h) results up to -1 scale height and then using the Chapman profile to peak. At 
night, the ionization at the lower altitudes not determined by the ionogram reduction technique more 
than offsets this effect. It is more accurate at night (and easier) to assume that the region is Chap- 
man rather than to integrate the profile obtained by the reduction. 

Figure 5 shows some variation in N mnx F 2 obtained by the different measurements. At the time of 
the rocket firing, the value ofN max F 2 at Fort Belvoir by the Chapman-fit method was 5.13 x 10 5 elec- 
trons/cm 3 , whereas the Wallops Island ionosonde gave 5.57 x 10 5 electrons/cm 3 . This indicates a 
horizontal gradient in the electron density, increasing to the east. The rocket firing was approxi- 
mately east of the ionosonde, and the measured value was 5.60 x 10 5 electrons/cm 3 . If the observed 
f 0 F 2 at Wallops Island is assumed to be correct, the horizontal gradient would have to be in the oppo- 
site direction east of the ionosonde. 


Another interesting observation is that the electron density profiles both day and night show a 
departure from the Chapman profile close to 225 km. Also, if the data from the daytime ionograms 
are used to calculate h max F 2 -H or, from the nighttime ionograms, h max F 2 - 2H , a value around 230 km 
is usually obtained. It would thus seem that the strong tendency to assume a Chapman distribution 
exists above about 225 km but not below it. 
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The ionogram data for July 13, 1960, was not included in Table 2. Although the day was not clas- 
sified as stormy, the F 2 region was undergoing rather rapid changes with time. None of the Fort Bel- 
voir ionograms near the time of the rocket firing agreed with the rocket data. The values obtained 
for H varied at times at a rate of 1 km per minute; in a few cases only three points were available 
within 1 scale height of the peak. Under such circumstances one cannot conclude that the profile was 
truly Chapman. It was noted, however, that the values of f Q F 2 computed on the assumption that it was 
Chapman gave values that were from 1 to 2 percent larger than the observed f Q F 2 . For these instances 
the recorded virtual heights were abnormally large, sometimes exceeding 600 km — presumably due 
to the fact that f Q F 2 was only slightly larger than f 0 F r This results in considerable absorption and a 
rapid rate of change in virtual height with frequency near f Q F 2 , which raises the question of whether 
it is possible to measure f 0 F 2 accurately from such an ionogram. Additional evidence is seen in the 
fact that during the period 20 minutes before firing time to 80 minutes after - during which ten iono- 
grams were taken — only one ionogram (20 minutes after firing) showed f Q F 2 at Fort Belvoir to be 
greater than the rocket measured value at Wallops Island, with only two (30 and 40 minutes after fir- 
ing) being nearly as large. 

If the ionogram reduction shows a Chapman distribution, one can obtain a reasonably accurate 
value for N b . Then, if measurements are made nearby of the total electron content by means of sat- 
ellites or two-frequency moon echoes, it is a simple matter to compute the value of the constant ex- 
ponential slope above h max F r The rocket values of N a + N b were assumed to be known, and the calcu- 
lation of these slopes agreed quite well in every case. It thus seems feasible to make measurements 
of the slope by this technique under quiet conditions, provided that it is not done near sunrise or 
sunset. 

The ratios N a /fa b that were obtained are quite interesting. Spring afternoon, summer morning, 
and summer evening values are all fairly close to 2.0. Garriott (Reference 13) also finds a similar 
value for autumn and winter days. On autumn and winter nights, however, he finds the ratio to be as 
high as 4. His two winter evening measurements made at 2044 local time are in very good agreement 
with the value of 3.2 reported here for that time on November 9. The model suggested in this paper 
shows that the ratio will increase if the value of Hnear h max decreases to lower values with no signif- 
icant change occurring in the constant logarithmic slope above h ma ^ As the ionosphere closely ap- 
proximates a Chapman region on the lower side at night, and as the ratio for a true Chapman region 
is 2.14 (Wright, Reference 1), one would expect that the ratio would be greater than this whenever the 
value of H is less than one-half the exponential slope, that is, about 75 km in the evening. 

Regarding disturbed ionospheric conditions, only a few remarks can be made. Berning* reported 
on a flight to 400 km at Wallops Island under disturbed conditions. The frequency f 0 F 2 was below 
normal, and h max F 2 was higher than normal. The F 2 region was very close to a Chapman distribution 
from -1/2 scale height to the peak of the flight, +1.3 scale heights. The value f 0 F t was nearly equal 
to f 0 F 2 , and under such conditions it is difficult to obtain an. ionogram reduction close to the peak. 
Also, the shape of the profile above +1.3 scale heights is not yet known in such cases. 


•See footnote, page I. 
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CONCLUSIONS 

The N(h) reduction of ionograms can be checked easily to determine whether the F 2 region fol- 
lows an a- Chapman distribution. If it does, the scale height H, the altitude F 2 , and electron con- 
tent N b can be determined with more consistent accuracy by a best-fit to a Chapman function than by 
using the parabolic technique. In addition, it is possible to predict accurately the electron density 
profile to 1.5 H above h max F 2 and to predict approximately the profile above this altitude, and there- 
fore to obtain a fair value for the electron content N a during the daytime and evening. If the total con- 
tent is measured by some method, the electron distribution may be predicted to nearly 1000 km by 
use of ionogram data and the suggested model, with the required logarithmic slope obtained from the 
total content measurement. 
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